Loop-Cluster Simulation of the t-J Model on the Honeycomb Lattice 
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Inspired by the lattice structure of the unhydrated variant of the superconducting material 
Nai:Co02-yH20 at a; = 5, we study the t-J model on a honeycomb lattice by using an efficient 
loop-cluster algorithm. The low-energy physics of the undoped system and of the single hole sector 
is described by a systematic low-energy effective field theory. The staggered magnetization per spin 
Ms = 0.2688(3), the spin stiffness ps — O.I02(2)J, the spin wave velocity c = 1.297(16) Ja, and the 
kinetic mass M' of a hole are obtained by fitting the numerical Monte Carlo data to the effective 
theory predictions. 

PACS numbers: 12.39.Fe, 75.10.Jm, 75.40.Mg, 75.50.Ee 



I. INTRODUCTION 

Since the discovery of high-temperature superconduc- 
tivity in cuprate materials [l|, the Hubbard and t-J mod- 
els have been of central importance in strongly corre- 
lated electron systems. However, due to their strong 
coupling, a systematic analytic treatment of these mod- 
els is currently not available. Similarly, a severe sign 
problem away from half-filling prevents us from under- 
standing these systems quantitatively by reliable Monte 
Carlo calculations. Despite these difficulties, much effort 
has been devoted to understanding the properties of t-J- 
type models on the square lattice. Although some con- 
troversial results have been obtained, various studies in- 
cluding exact diagonalization |2,|3|, series expansion [J|, 
and Monte Carlo simulations [5|, |6| enable us to under- 
stand the hole-dynamics quantitatively, at least to some 
extent. In particular, these studies all obtained minima 
of the single-hole dispersion relation at lattice momenta 
(±-^,±-^) in the Brillouin zone of the square lattice 
which is in agreement with experimental results [7|, |8|, |9[ . 

A reliable and order-by-order exact way to investi- 
gate the low-energy physics of lightly doped antiferro- 
magncts is provided by a systematic low-energy effec- 
tive field theory. The physics of the undoped systems 
is quantitatively described by magnon chiral perturba- 
tion theory [l^, [HI, Ej E, [IJIj while the interactions 
of magnons and holes are described by a low-e nerg y ef- 
fective theory for hole-doped antiferromagnets [I^, |T^ . 
Predictions of the effective theory only depend on a small 
number of low-energy constants which can be determined 
from cither experiments or Monte Carlo data. Thus, the 
use of low-energy effective theories together with reliable 
Monte Carlo simulations provides an unbiased approach 
to studying the low-energy physics of these systems. In 
particular, using the loop-cluster algorithm [17[, the low- 
energy parameters of the spin 1/2 Hciscnberg model have 
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been determined with very high precision |18l . [lOf. In- 
deed, thanks to the combination of very efficient Monte 
Carlo simulations with the systematic low-energy effec- 
tive field theory, undoped antiferromagnets on the square 
lattice like La2Cu04 and Sr2Cu02Cl2 are among the 
quantitatively best understood condensed matter sys- 
tems. 

In addition to the cuprates, another superconducting 
material, Na2,Co02-yH20, has drawn a lot of attention 
both theoretically and experimentally. Unfortunately, 
due to the fact that the underlying lattice geometry of 
the spin 1/2 cobalt sites in these materials is triangular 
— which leads to strong geometric frustration — a first 
principles Monte Carlo study is impossible in practice. 
Nevertheless, the spin- and charge-ordering tendencies 
observed and studied in [20, UM, ^M ™^y suggest that at 
filling X ~ ^, the unhydrated parent compound Naj;Co02 
can be described by the t-J model on a half-filled hon- 
eycomb lattice which allows one to simulate the system 
efficiently with the loop-cluster algorithm. 

Another system on the honeycomb lattice that has 
been investigated with great vigor is graphene — a sin- 
gle sheet of graphite (see [23| for a detailed review) . As 
a consequence of the geometry of the honeycomb lattice, 
the low-energy excitations of graphene are massless Dirac 
fermions. If some variants of graphene exists at stronger 
coupling, one eventually expects a phase transition sepa- 
rating graphene's unbroken phase from a strong coupling 
antiferromagnetic phase in which the SU{2)s spin sym- 
metry is spontaneously broken to U{l)s- The low-energy 
effective theory of the unbroken phase and of the critical 
point has been constructed in [24| . 

Motivated by possible applications to Naa;Co02, we in- 
vestigate the spin 1/2 Heisenbcrg model as well as the t-J 
model on the honeycomb lattice by using the quantum 
Monte Carlo method. Just as in the square lattice case, 
the long-distance physics of these models is described 
quantitatively by a systematic low-energy effective field 
theory. At low energies, the Heisenberg model on a bi- 
partite lattice is described by magnon chiral perturbation 
theory, and, accordingly, the t-J model is described by a 
low-energy effective field theory for magnons and holes. 
Based on the same method that has been used in the 



square lattice case, we have constructed the leading-order 
terms in the action of a systematic low-energy effective 
field theory for magnons and holes on the honeycomb 
lattice. In this paper, we determine the corresponding 
leading-order low-energy constants, namely the staggered 
magnetization per spin Ms, the spin stiffness ps, the spin 
wave velocity c, and the kinetic mass M' of a doped hole 
by fitting the Monte Carlo data to the effective field the- 
ory predictions. 

The rest of this paper is organized as follows. In sec- 
tion |lll we introduce the relevant microscopic models as 
well as corresponding observables. Section|TTT]reviews the 
low-energy effective theory for magnons, and section IIVI 
is devoted to the Monte Carlo determination of the corre- 
sponding low-energy parameters. The single-hole physics 
is investigated in section |Vl and the effective theory for 
holes and magnons is discussed in section IVIl Finally, 
section IVTIl contains our conclusions. 



II. MICROSCOPIC MODELS AND 
CORRESPONDING OBSERVABLES 

In this section we introduce the Hamiltonians of the 
microscopic t-J model and the Heisenberg model as well 
as some relevant observables. The t-J model is defined 
by the Hamilton operator 



H = P{- t^iclcy + do 



{xy) (xy) 



Here c|, and Cx are fcrmion creation and annihilation op- 
erators at a site x with 



Cxi 



(2) 



whose components obey standard anticommutation rela- 
tions. In terms of the Pauli matrices a the local spin 
operator at a site x is given by 



^x 



J 



(3) 



The projection operator P restricts the Hilbcrt space by 
eliminating doubly occupied sites. Hence the t-J model 
allows empty or singly occupied sites only. The hop- 
ping of fermions is controlled by the parameter i, while 
J > is the antiferromagnetic exchange coupling be- 
tween neighboring spins. At half-filling, the t-J model 
reduces to the Heisenberg model with the Hamiltonian 



H 



{xy) 



(4) 



The honeycomb lattice with periodic spatial boundary 
conditions implemented in our simulations is depicted in 
figure 1. The dashed rectangle in figure 1, which contains 
4 spins, is the elementary cell for building a periodic hon- 
eycomb lattice covering a rectangular area. For instance. 



the honeycomb lattice shown in figure 1 contains 3x3 
elementary cells. The lattice spacing a is the distance 
between two neighboring sites. The honeycomb lattice is 
not a Bravais lattice. Instead it consists of two triangu- 
lar Bravais sub-lattices A and B (depicted by solid and 
open circles in figure 1). As a consequence, the momen- 
tum space of the honeycomb lattice is a doubly-covered 
Brillouin zone of the two triangular sub-lattices (depicted 
in figure 2). 

A physical quantity of central interest is the staggered 
susceptibility which is given by 



1 



Xs 



P 



L1L2 Jo 

1 /-/^ 



dt (M!{0)M!it)) 



rfiiTr[Mf(0)M3(i)exp(-/3iJ)].(5) 

^l-'^2 Jo '^ 

Here /3 is the inverse temperature, Li and L2 are the 
spatial box sizes in the xi- and a;2-direction, respectively, 
and 



Z = Trexp(-/3iJ) 



(6) 



is the partition function. The staggered magnetization 
order parameter Mg is defined by 



Ms = 5](-lf ^. 



(7) 



Here (-1)^ = 1 on the A- and (-1)^ == -1 on the B- 
sub-lattice, respectively. Another relevant quantity is the 
uniform susceptibility which is given by 

Xu = -r-^ I dt (A/3(0)Af3(t)) 



L^L 



1-^2 Jo 



1 



L1L2 



'^ 1 
dt -Tr[M3(0)M3(i) exp(-/3i/)]. (8) 



Here 



M 



Y.sx 



(9) 



is the uniform magnetization. Both Xs Eind Xu can be 
measured very efficiently with the loop-cluster algorithm 
using improved estimators [l8| . In particular, in the 
multi-cluster version of the algorithm the staggered sus- 
ceptibility is given in terms of the cluster sizes \C\ (which 
have the dimension of time), i.e. 



1 



4/3L2 



Eici^ 



Similarly, the uniform susceptibility 



Xn 



4L2 



{W?) 




Y.^t{Cf 



(10) 



(11) 



is given in terms of the temporal winding number Wt = 
5^^ Wt{C) which is the sum of winding numbers Wt{C) 
of the loop-clusters C around the Euclidean time direc- 
tion. Similarly, the spatial winding numbers are defined 
by Wi = J2C ^*(<^) with i e {1, 2}. 




FIG. 1: The periodic honeycomb lattice consisting of two tri- 
angular sub-lattices A and B, which are depicted by solid and 
open circles, respectively. The dashed rectangle is an elemen- 
tary cell for building a periodic honeycomb lattice covering a 
rectangular area. 




FIG. 2: The momentum space of a honeycomb lattice, which 
is a doubly- covered Brillouin zone dual to the two triangular 
sub-lattices A and B. 



III. LOW-ENERGY EFFECTIVE THEORY FOR 
MAGNONS 

Due to the spontaneous breaking of the SU{2)s 
spin symmetry down to its C/(l)s subgroup, the low- 
energy physics of antifcrromagncts is governed by two 
massless Goldstonc bosons, the antiferromagnetic spin 
waves or magnons. The description of the low-energy 
magnon physics by an effective theory was pioneered by 
Chakravarty, Halperin, and Nelson in [10|. In analogy to 
chiral perturbation theory for the pseudo-Goldstone pi- 
ous in QCD, a systematic low- energy effective field the- 
ory for magnons was developed in [ll|, [i3, ll3, 13 • The 



staggered magnetization of an antiferromagnet can be 
described by a unit-vector field e{x) in the coset space 

5f/(2),/C/(l), = 52, i.e. 



e{x) = (ei(x),e2(x),e3(a;)), e{xf 



1. 



(12) 



Here x = {xi, X2, t) denotes a point in (2+l)-dimensional 
space-time. To leading order, the Euclidean magnon low- 
energy effective action takes the form 



S[e] 



d X dt 



Ps 



die- die 



1 



dte- dte 



(13) 



where the index i € {1,2} labels the two spatial direc- 
tions and t refers to the Euclidean time-direction. The 
parameter ps is the spin stiffness and c is the spin wave 
velocity. At low energies the antiferromagnet has a rela- 
tivistic spectrum. 

Using the above Euclidean action, detailed calculations 
of a variety of physical quantities including the next-to- 
next-to-leading order contributions have been carried out 
in [2^. Here we only quote the results that are relevant 
for our study, namely the finite-temperature and finite- 
volume effects of the staggered and uniform susceptibil- 
ities, as well as results on the rotor spectrum of the an- 
tiferromagnet in a finite volume. The aspect ratio of a 
spatially quadratic space-time box with Li ~ L2 ~ L 
is characterized by I = {(3c/ L)^/^ , with which one dis- 
tinguishes cubical space-time volumes with fSc k. L from 
cylindrical ones with f3c ^ L. In the cubical regime the 
volume- and temperature-dependence of the staggered 
susceptibility is given by 



MiL^P 



1 + 2- 



7/3i(0 



PsLl 



PsLl 

[/3i(0' + 3/32(0] +0 




,(14) 



where A4s is the staggered magnetization density. The 
uniform susceptibility takes the form 



ll{'-5;fll*<" 



1 

3 
O 



PsLl 
1 

Z3 



/32(0 



-(3,ilf^6iil) 



(15) 



The functions (}i{l), Pi{l), and V'(0, which only depend 
on Z, are shape coefficients of the space-time box defined 
in [25|. In the very low temperature limit, one enters the 
cylindrical regime of space-time volumes with /3c ^ L. 
In this case, the staggered magnitization vector M^ acts 
as a quantum rotor and, correspondingly, the low-energy 
end of the spectrum takes the form 



Es 



S{S + 1) 
26 



(16) 



Here S G {0, 1, 2, ...} is the spin and is the moment of 
inertia of the quantum rotor which is given by [25| 



6 = 



PsL^ 



3.900265c ^ / 1 



(17) 



The partition function of the (25+l)-fold degenerate ro- 
tor spectrum is given by 



Z = J2{'2S + l)cxp{-f3Es). 



(18) 



s=o 



The probabihty distribution of the uniform magnetiza- 
tion M^ = S^ is then given by 



<M') = ^ E exp(-/3i?s). 



(19) 
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IV. DETERMINATION OF THE LOW-ENERGY 
PARAMETERS OF THE UNDOPED SYSTEM 

In order to determine the low-energy constants Ms, 
Ps, and c, we have performed numerical simulations 
of the Heisenberg model on the honeycomb lattice 
with up to 4680 spins in the cubical and cylindrical 
regimes. The cubical regime is determined by the condi- 
tion iEcWiiCf) « (Ec^2(C)2) « (EcWtiCr) 
(which implies (3c ~ L). The chiral perturbation theory 
predictions for Xs and Xu in Eqs. (jl4|) and (|15[) are derived 
for a (2-|-l)-dimensional box with equal extent in the two 
spatial directions (which we refer to as a square-shaped 
area) . Since it is not possible to consider the honeycomb 
geometry on an exactly square-shaped area, our simula- 
tions are done on almost square-shaped rectangles. To 
be more precise, the lattices used in our simulations de- 
viate from a perfect square-shaped area by less than 0.4 
percent. We have performed an interpolation on some 
of our data to the exactly square-shaped area and find 
agreement between the fits of the interpolated data and 
the raw data. The inclusion of 0{1/L^) corrections in 
the fits leads to consistent results as well. Instead of 
considering the staggered magnetization density Ais of 
Eg. p^ . we choose to quote the staggered magnetization 
per spin Ms, which is related to A4s by 



Ms 



3^3 



Msa^ 



(20) 



TABLE I: Some numerical data for the staggered susceptibility 
Xs and the temporal winding number squared {Wf ) obtained 
with the loop-cluster algorithm. Ni and N2 count the number 
of copies of elementary rectangles in the 1- and 2-direction 
and Nspin = 4N1N2 is the corresponding number of spins. 



the permille level. The value of c obtained here is con- 
sistent with the one of a spinwave expansion study [26|. 
The above value oi Ms is larger than the one of a pre- 
vious spinwave expansion 12711 but consistent with that 
of a series expansion study [28| (within the comparably 
large 4 percent error of that study). It is only slightly 
larger than the value obtained in a previous Monte Carlo 
calculation Ms = 0.2677(6) ^. We want to point out 
that our results are obtained by fitting more than 80 nu- 
merical data points to two analytic predictions with only 
3 unknown parameters. If Ms is fixed to 0.2677, the 
quality of our fit downgrades to x^/d.o.f. ~ 3.0. The re- 
duction of Ms = 0.2688(3) and ps = 0.102(2) J on the 
honeycomb lattice compared to those on the square lat- 
tice {Ms = 0.3074(4), Ps = 0.186(4) J [ll,[Il) indicates 
larger quantum fluctuations on the honeycomb lattice. 
This is expected since the coordination number of the 
honeycomb lattice is smaller than the one of the square 
lattice. 



Some numerical data from our simulations arc listed in 
table 1. By fitting Xs and Xm simultaneously to Eq. fTil) 
and Eq. (fT5|) . we find 

Ms = 0.2688(3), Ps = 0.102(2) J, c = 1.297(16) Ja (21) 

with xVd-o.f. « 1.05 (see figures [3] and E]) . The low- 
energy constants ps and c are determined with high ac- 
curacy (at the percent level). The error oi Ms is even at 



Having determined the values of the low-energy param- 
eters Ms, Ps, and c from the cubical space-time volume 
regime, we can test the effective theory in the cylindrical 
regime. Figure [5] shows a comparison of the effective the- 
ory prediction for the probability distribution p{M^) of 
Eg.pg]) with Monte Carlo data. The observed excellent 
agreement — which does not involve any adjustable pa- 
rameters — confirms the quantitative correctness of the 
effective theory. 




FIG. 3: Fit of the finite-size and finite-temperature effects of 
the staggered susceptibility Xa to results of the effective theory 
in the cubical regime. 
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FIG. 4: Fit of the finite-size and finite-temperature effects of 
the temporal winding number squared (Wt) to the results of 
the effective theory in the cubical regime. 



V. SINGLE HOLE DISPERSION RELATION 

The physics of a single hole on the honeycomb lat- 
tice was studied theoretically in [30| using exact diago- 
nalization, series expansion, and self-consistent Born ap- 
proximation. Here we use first principle Monte Carlo 
simulations to quantitatively investigate the single-hole 
dispersion relation and quasiparticlc weight. To achieve 
this goal, we have implemented a technique similar to 
the ones used in [^, Q to simulate the one-hole sector of 
the t-J model. To calculate the fcrmion two-point func- 
tions in momentum space, one should keep in mind that 
the Brillouin zone of the honeycomb lattice is doubly 
covered. The two covers are dual to the two triangular 
Bravais sub-lattices A and B. Therefore, one needs to 
distinguish the correlators between AA, AB, BA, and 
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FIG. 5: Comparison of the effective theory prediction for the 
probability distribution p{M^) of Ea. Uy\) with Monte Carlo 
data on an Ni X N2 = 11 x 19 honeycovfib lattice with Nsptn = 
836 at pj = 60. The open circles are the Monte Carlo data 
while the histogram is the effective theory prediction. 



BB sub-lattices. The correlation function between AA 
sub-lattices with momentum k takes the form 



G^^(fc,t) = I E Tr[ct(0)c,(i)exp(-/3i/)] 

x,y£A 

X exp(— iA:(a; — y)) 



Y, Zn{k) cxp (-(S„(fc) - Ea)t) , (22) 



where Eq is the ground state energy of the half-filled 
system, and 



Zn{k) = |(0| E Cxex.v{ikx) 



(23) 



The factor Zi{k) is known as the quasiparticlc weight. 
In deriving Eg. ip^ . we have inserted a complete set of 
energy eigenstates 1 = ^„ |"')("'| in the single-hole sector 
and taken the limit /? — s- 00 in the final step. The fcrmion 
energy 



Eh{k) ^ E^{k) - E^ 



(24) 



corresponding to the momentum k can be extracted by 
fitting the data to a single- or a double-exponential. The 
correlation function between AA sub-lattices with mo- 
mentum A: = (f^, J^J- ) depicted in figure [H] is obtained 
on a honeycomb lattice with 3456 spins and J /t = 2.0. 
A single-exponential fit yields Eh{k) — 0.207(9)t while a 
double-exponential fit results in Eh{k) = 0.201(5)i. The 
two fits yield consistent results. In the same way, we 
determine the one-hole dispersion relation from the AA 



FIG. 6: Correlation function G^^{k,t) between A A sub- 
lattices with Fourier momentum k = ('^, „'^%- ). The bottom 
line is the result of a single-exponential fit while the top line 
is obtained from a double- exponential fit. 



lattice has been constructed using the known location of 
hole or electron pockets. This powerful method was used 
to systematically construct the effective theory for t-J- 
type models on the square lattice in [l^, [l^, [Sa] ■ The ef- 
fective theories were used to investigate the one-magnon 
exchange potential and the resulting bound states be- 
tween two holes and two electrons as well as the possible 
existence of spiral phases [ifl, |3g, [sj ■ Using the informa- 
tion about the location of the pockets and based on the 
symmetry properties of the underlying microscopic the- 
ory, we have constructed a systematic low-energy effec- 
tive theory for the t-J model on the honeycomb lattice. 
The details of the construction of the effective theory 
will be described in a forthcoming publication. Here we 
briefly sketch the principles behind this construction and 
present the terms in the effective Lagrangian that arc rel- 
evant to our present study. In the effective theory, the 
holes which reside in momentum space pockets centered 
at 



fc" = (0, 



47r 
3\/3a 



47r 



^•"'-^"■-JA^ 



(25) 




FIG. 7: The dispersion relation Eh{k)/t for a single hole in 
an antiferromagnet on the honeycomb lattice. 



correlator for all momenta k. The single-hole dispersion 
relation in figure [7] is obtained with the same parameters 
as in figure [5] The figure shows that the hole pockets 
are located at (±|f ,±^^) and (0,±^i|^) in the Bril- 
louin zone. The position of the hole pockets agrees with 
the position of the Dirac cones obtained from the free 
fermion theory on the honeycomb lattice which is rele- 
vant for graphene. 



are represented by Grassmann fields ipl{x). Here the 
"flavor" index f = a,f3 characterizes the corresponding 
hole pocket and the index s = ± denotes spin parallel (-I-) 
or anti-parallel (— ) to the local staggered magnetization. 
The magnons are coupled to the holes through a nonlin- 
ear realization of the spontaneously broken SU{2)s sym- 
metry. The global SU{2)s symmetry then manifests it- 
self as a local U{l)s symmetry in the unbroken subgroup. 
This construction leads to an Abelian "gauge" field w|J (x) 
and to two vector fields v^ {x) which are "charged" under 
U{l)s spin transformations. The coupling of magnons 
and holes is realized through vf^{x) and v^{x). These 
fields have a well-defined transformation behavior under 
the symmetries which the effective theory inherits from 
the underlying microscopic models. Based on symme- 
try considerations, we have constructed the leading or- 
der terms of the effective Lagrangian for magnons and 
holes on the honeycomb lattice. In this paper we only 
list those terms that are relevant for the propagation of 
a single hole, i.e. 






(26) 

Here M is the rest mass and M' is the kinetic mass of a 
hole, while I?^ is a covariant derivative given by 



VI. EFFECTIVE FIELD THEORY FOR HOLES 

AND MAGNONS 

Motivated by baryon chiral perturbation theory for 
QCD [3l|, [33, [33, [SJ, [3a], based on symmetry considera- 
tions, a universal effective theory for magnons and charge 
carriers in lightly doped antiferromagnets on the square 



D^i;i{x) = [^,±^vl{x)]i;i{x). 



(27) 



Eq. (p6|) yields circular hole pockets for small momenta 
which is indeed confirmed in figure [H] The low-energy 
constant M' in Eq. ((26|) is obtained from the curvature of 
the dispersion Eh{k) near a minimum. For example, on 
a honeycomb lattice with 3456 spins and J /t = 2.0, we 
find A/':=4.1(l)/(ta2). 
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TABLE II: Kinetic mass M' as well as the bandwidth A for 
some values of J/t. 
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FIG. 8: Circular hole pockets on the honeycomb lattice. The 
dot corresponds to the point Gamma and the stars mark the 
centers of the hole pockets (corresponding to the point W and 
its symmetry partners). The parameters are the same as in 
figure^ 



In figures [9] and [TOl we have plotted the single-hole 
dispersion as well as the quasiparticle residue Zi{k) over 
the irreducible wedge T-K-W-T of the Brillouin zone for 
J/t = 1.0. The resulting bandwidth 



A ^ Eh{r) - Eh{W) 



(28) 



is in qualitative agreement with exact diagonalization 
and series expansion in [30|. While exact diagonaliza- 
tion of small systems may suffer from finite size effects, 
and series expansions may not converge in all regions of 
parameter space, the Monte Carlo data obtained with the 
efficient loop-cluster algorithm do not suffer from system- 
atic uncertainties. In table 2 we list the kinetic mass M' 
as well as the bandwidth A for a few values of J/t. 



VII. CONCLUSIONS 

We have studied antiferromagnetism on the honey- 
comb lattice by first principle Monte Carlo calculations. 
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FIG. 9: Dispersion relation Eh{k)/t of a single hole for J/t = 
1.0 along the irreducible wedge T-W -K-V m the first Brillouin 
zone (see figure 2). 




FIG. 10: The quasiparticle weight Zi{k) of a single hole for 
J/t = 1.0 along the irreducible wedge F-W-K-F in the first 
Brillouin zone (see figure 2). 



In particular, we have fitted more than hundred Monte 
Carlo data obtained at rather different volumes and tem- 
peratures with just four low-energy parameters — Ms, 
Ps, c, and M' — of the effective theory for the t-J model. 
These parameters have been determined with percent 
and sometimes even with permille accuracy. This should 
demonstrate convincingly that the systematic low-energy 
effective field theory yields quantitatively correct results 
for the physics of magnons and holes. Hence, the effec- 
tive theory allows us to perform unbiased investigations 
of the low-energy physics of the system. The construction 
of the effective theory for lightly doped antiferromagnets, 
as well as a systematic investigation of the one-magnon 
exchange potential and the resulting bound states be- 
tween two holes as well as the possible existence of spiral 
phases of lightly doped antiferromagnets on the honey- 
comb lattice will be presented in subsequent studies. 
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